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ABSTRACT 

v  \ 

He  consider  boundary  value  problems  for  the  biharmonic  equation  in  the  open 
rectangle  x  >  0,  -1  <  y  <  1 ,  with  homogeneous  boundary  conditions  on  the  free 
edges  y  *  fl ,  and  data  on  the  end  x  »  0  of  a  type  arising  both  in  elasticity  and 
in  Stokes  flow  of  a  viscous  fluid,  in  which  either  two  stresses  or  two  displacements 
are  prescribed.  For  such  " non-canon ical“  data,  coefficients  in  the  eigenfunction 
expansion  can  be  found  only  from  the  solution  of  infinite  sets  of  linear  equations, 
for  which  a  variety  of  methods  of  formulation  have  been  proposed. 


A  drawback  of  existing  methods  has  been  that  the  resulting  equations  are 
unstable  with  respect  to  the  order  of  truncation.  It  is  clear  from  an  examination 
of  the  spectrum  of  a  typical  matrix  that  ill-conditioning  is  to  be  expected. 

However,  a  search  among  a  wider  class  of  possible  trial  functions  than  hitherto  for 
use  in  a  Galerkin  method  based  on  the  actual  eigenfunctions  has  led  to  the  choice  of 
a  unique  set,  here  termed  optimal  weighting  functions,  for  which  the  resulting 
infinite  matrix  is  diagonally-dominated.  This  ensures  the  existence  of  an  inverse, 
which  can  be  approximated  by  solving  a  finite  subset  of  the  equations, 


Computations  for  a  number  of  representative  cases,  presented  in  full  in  an 
internal  report  (Spence  1978)  are  summarized  here,  with  emphasis  on  the  rates  of 
decay  of  the  coefficients  {c  }  in  the  eigenfunction  expansion.  Knowledge  of  these 
decay  rates  is  essential  for  a  discussion  of  convergence,  parallel  to  .nat  given  by 
Joseph  ( 1977a, b)  and  his  co-workers  for  canonical  problems. 


Asymptotic  estimates  of  the  decay  rates  have  also  been  obtained  by  use  of  the 
solution  of  the  biharmonic  equation  in  a  quarter  plane.  It  is  found  that  (i)  for 
smooth  continuous  data  satisfying  compatibility  conditions  at  the  corners,  the  decay 
rates  guarantee  pointwise  convergence.  Also  examined  are  (ii)  cases  of  data 
violating  compatibility  (iii)  discontinuous  data  and  (iv)  discontinuities  in 
derivatives  of  the  data.  In  these  cases  sharp  estimates  of  convergence  rates  are 
obtained,  which  guarantee  that  integrals  of  the  series  converge  to  integrals  of  the 
data.  The  computations  show  striking  confirmation  of  the  theoretical  estimates. 
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A  CLASS  OF  BIHARMONIC  END-STRIP  PROBLEMS  ARISING 
IN  ELASTICITY  AND  STOKES  FLOW 


* 

D.  A.  Spence 


1.  Introduction 

The  solution  of  the  bihsmonic  equation 


A2  Y  -  0 


(1.1) 


in  a  semi-infinite  strip  fl:  {  *  £  0,  |y|£  1}  with  homogeneous  boundary 
conditions  on  the  edges  y  ■  +  1  can  be  expressed  formally  as 

*  (*»y)  -  l  cn  i|»n(y)e"Xn*  (1.2) 

where  $  (y)  =  (y,X  )  are  eigenfunctions  associated  with  the  names 
n  n 

of  Papkovitch  (1940)  and  Fad'le  (1941),  Xq  are  eigenvalues,  and  {cq}  are 

coefficients  that  depend  on  the  data  on.  the  edge  T:  {x  -  0,  |y|£  l). 

In  the  case  of  solutions  even  in  y,  and  such  that  f  +  0  as  x  +  »,  with  edge 
conditions 


9  ■  0,  T  ■  0  on  y  -  +  1,  (1.3) 

the  eigenvalues  are  the  zeros  with  positive  real  parts'  of  the  function 


C(X)  5  X  ♦  sin  X  cos  X  (1.4) 

and  the  eigenfunctions**  may  be  taken  as 

)p(y,X)  “  (X  cos2X)  ^ (sinX  cosXy  -  y  cosX  sin  Xy) 

(1.5) 

=  (X  cos2X)’1  $  (y,X)  say 


*Mathematics  Department,  Imperial  College,  London  SW7  2BZ,  England. 

**These  satisfy  the  equation  (D2+A2)2*  -  0  (D  h  )  with  boundary  conditions 
<»(1)  ■  D^(l)  ■  0,  D$(0)  -  D3ip (0)  -  0.  They  are  alsoydiscussed  by  Lure  (1964) 
Buchwald  and  Doran  (1964)  and  many  other  authors  in  addition  to  those  cited. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


The  present  paper  is  concerned  only  with  data  that  can  be  expanded  in 
terms  of  these  even  eigenfunctions.  However  there  is  also  a  set  of 
eigenfunctions  that  are  odd  with  respect  to  y.  the  eigenvalues  in  that  case 
being  seros  of 


C  (X)  =  X  -  sinX  cosX 


and  the  eigenfunctions  multiples  of 


(1.6) 


♦  =  cosX  sinXy  -  y  sinX  cos  Xy 


(1.7) 


In  general,  an  expansion  such  as  (1.2)  would  involve  eigenfunctions  of  both 


sets. 


Footnote:  Different  authors  have  used  different  normalising  factors 

from  the  present  (Xcos2X)  *  when  defining  coefficients  in  the  expansion 
(1.2  )  in  terms  of  the  basic  eigenfunctions  $.  Smith  (1952),  who  is  followed 
hy  Joseph  (1977  et  seq)  writes  the  expansion  as  F  ■  e~\iX  ♦(y.X^), 

while  Gregory  (1980  a.b)  writes  F  -  l  (-2)c e*\»x  $(y,|X  )  (so  that  his 

n  n 

Xn'*  *r*  twice  those  defined  by  (1.4),  i.e.  they  are  the  roots  of  X  ♦  sinX  »  0) , 

Thus  the  coefficients  c  obtained  in  the  present  paper  are  cos2X  times  those 

n  o 

of  Joseph,  and  -  |Xn  cos2Xn  times  those  of  Gregory.  For  n  »  1,  these  ratios 
are  0(n)  and  0(n2)respectively. 

o  f\  Th*  reason  for  the  present  choice  of  normalizing  factor  is  that  it 
S  9  •  <v  I 

*1  /'!*•<*•  to  the  particularly  simple  biorthogonality  relations  (2.9  ),  (2.10). 


'  -*r~7  ~  r 


I  -■)  • 


1  I  f  w 


C(X)  has  a  siaple  zero  at  X  -  0,  and  in  addition  a  conjugate  pair  of 

zeros  X  ,  X  in  each  interval  (n-i)ir  <  ReX  <  (n-l/4)ir.  -X  ,  -X  are 
n  n  n  n 

also  zeros.  Since  we  are  concerned  in  this  paper  only  with  those 
zeros  with  positive  real  parts,  we  shall  adopt  the  numbering  convention 
that  X^,  X?  .  are  the  zeros  in  the  first  quadrant,  and 


n  ■  1 ,  2,  3,  ....  . 


(1.8) 


The  sum  (1.2)  is  then  to  be  understood  as 


l  cn  4>n(y)  e  n 


<n#» 


_  ,  _  c  .  -X  x 

=  2  Re  I  cn  V  n  * 

n  ■  1 


the  zero  eigenfunction  being  excluded. 

The  asymptotic  location  of  the  roots  is  given  by 


X  «  (n  -  l/4)ir  ♦  In  4mr  +  o  (1) 
n  2 


This  expression  shows  that  as  n  -*•  °» 


while 


Likewise 


|  X  J  'v.  nx,  j  cosX  j,  |sinX  |  (nir)^ 


tan  X  -  i  (1-  1/4X  2)1  +  i. 
n  n 


tan  X  -*•  -  i 


(1.9) 


(1.10) 


(1.11) 


The  first  10  eigenvalues,  together  with  a  much  closer  asymptotic  expression, 
are  listed  in  Table  1. 

The  zero  eigenfunction  satisfying  the  edge  conditions  in  the  form 
9  ■  Vy  ■  0  is  identically  zero.  However  if  the  edge  conditions  are 
taken  as  •  0,  H^y  -  0,  which  correspond  in  elastostatics  to 
vanishing  tractions,  a  term  y2  corresponding  to  X^  could  be  included 
in  (1.9  ). 


Table  1.  First  10  stroi  of  2A  ♦  sin  2A 


2.10619617 

5.35626888 

8.53668213 

11.69917774 

14.85406017 

18.00493240 

21.15341377 

24.30034256 

27.44620323 

30.59129524 


1.12536430 

1.55157435 

1.77554369 

1.92940450 

2.04685259 

2.14189076 

2.22172284 

2.29055238 

2.35104823 

2.40501261 


Tha  asynptotic  expression  £  C  -  ln  Cn  gives  6  figure  agreement  with 

'n 

Ba  An  at  n  -  4.  and  {  ln  5n  +  [<ln  C  )2-21n  C-  l]/t  2  agrees  with  Is>  \ 


to  4  x  10  at  n  -  9  (cn  -  (4n  -  l)ir). 


■  *.vA*sai 


2  Boundary  value  problems 


The  determination  o£  the  coefficients  {c  }  in  the  expansion 

n 

(  1.2  )  requires  knowledge  of  two  sets  of  boundary  data  on  the  edge 
T:  x  ■  0,  |y|<,l.  For  this  purpose  following  a  numerical  classification 

scheaie  of  the  type  introduced  by  Johnson  &  Little  to  treat  a  range  of  mixed 
boundary  value  problems  we  may  define  the  following  quantities: 


*y 


(1)(y). 

^yy(O.y)  -  f(2)(y) 

<5><y>. 

P(0,y)  -  f(4)(y) 

Hare  P  (x,y)  is  the  Laplacian  A¥  and  Q  the  harmonic  conjugate  function 
defined  so  that  P  ♦  iQ  is  analytic  in  x  ♦  iy,  with  Q(0,0)  -0.  An 
alternative  notation  used  by  Smith  (1952)  and  followed  by  Joseph  Q.977  ) 
and  Gregory  (1980)  is  to  write 


*xx  <°»y>  -  f(y>»  'ryy(0,y>  “ 8(y) 


To  this  it  would  be  convenient  to  add 


*xy  *  h*  p  ■  P»  Q  "  9  on  x  -  0 


Then  f  S  f<4)  -  f<2\  g-f<2\  h-f(1>  p-f<4\  q-f<3>. 


(2.1) 


(2.2) 


(2.3) 

(2.4) 


-5- 


Fomally, 


the  derived  functions  possess  expensions  of  the  form 


and  with  x  *  0,  following  the  notation  of  the  author's  earlier  paper, 
the  expansions  can  be  written 


H 

-*■  \ 

1  '°l)\ 

f  (2) 

r  1 

1 

A2) 

♦n 

■  Ic 

o' 

K 

III 

f  O) 

u  n 

2sinAy/cosA 

f(^>  1 

p2cosAy/cosA  / 

X“Xn 

l  <‘7 

Alternatively  one  can  write  expressions  for  f,  g  and  h  as 


t*  \ 


\"  / 


h. 


while 


J2(c  /cosA  ) 
n  n 


-cosA 

sinA 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


aiffiiihiinnm 


-6- 


Canonical  problems 


The  determination  cf  the  constants  {c  }  in  the  expansion  (  1.2)  requires 

n. 

knowledge  of  two  of  the  functions  f^  (or,  more  generally,  of  two  linear 

combinations  of  the  f^).  In  general  the  {cq}  must  then  be  found  from  an 

infinite  set  of  linear  equations,  but  there  are  two  problems  for  which  they  can 

be  found  explicitly,  by  use  of  certain  biorthogonality  relationships  among  the 

},  namely 
n 


1  <V3)  *(1>  ♦  *(1V3)  >  -  6 

*  Ym  n  m  n  /  mn 

and  J  <><4>  *<2>  +  (*(2>V4)>*<4)>-  6 

2  \  m  Tn  m  Ym  'Tn  /  mn 

Here  \*  /  denotes  (•)  dy. 

Thus  (1)  for  the  problem  in  which  the  prescribed  data  is  the  boundary  values 

of  f^  and  f^3\  the  scalar  product  of  (2.6  )  with  ($^2\  0,  ,  0)  gives 

m  m 


(2.9 

(2.10 


c.  - 1  <♦«>*<»  ♦  ♦«>,<»> 

as  noted  in  Appendix  A, 

In  the  context  of  plane  elasticity/  with  Y  as  the  Airy  stress  function, 
corresponds  to  a  case  in  which  the  shear  stress  -  Y  and  the  normal 
displacement  u  are  prescribed  over  T. 

(2)  Likewise  the  product  with  (0,$^4\  0,  4^  -  4^)  gives 

d  in  m 


(2.11) 


this 


(2.12) 


In  the  same  way,  this  corresponds  to  the  normal  stress  Y  and  shear 

yy 

displacement  v  being  prescribed  over  F. 


In  Um  of  the  classificaticn  (2.7-8),  the  last  two  expressions 

are 


en  -  (1/2) 

\  (2sinXy/cosX)  h  - 

II 

(2.13) 

and 

cn  -  (1/2) 

^(-2cosXy/cosX)  8  -  (X2*)p)> 

n 

(2.14a) 

-  (1/2) 

<  -X’*f  *  ♦“&* 

Q 

(2.14b) 

respectively. 


* 

1 


3  Mixed  problen* 


Three  non-canonical  problems  of  fundamental  importance  can  be 
identified*  as  noted  in  Appendix  A. 

1.  The  elastic  end-stress  problem:  ^Xy*  Vyy  prescribed  on  T 

2.  The  elastic  end-displacement  problem:  ¥  ♦  (l-v)Q,  T  -(l-v)P, 

*jr  // 

(v  ■  Poisson's  ratio)  prescribed. 

3.  The  end-stress  problem  in  Stokes  flow  :  fxy+  JQ,  Tyy-|P 
prescribed  on  T 

All  three  of  these  problems  fall  into  the  same  framework  if  the  data 
is  prescribed  in  the  form  of  values  of  the  two  functions 

g<»  -  f<»  ♦  a  f<3>,  g<2>  -  f<2>  -  af(4> 

with  a  -  0,  1  -  v,  $  in  cases  1,  2,  3  respectively. 

For  these  problems,  the  biorthogonality  relations  cannot  be  used  to 
provide  explicit  expressions  for  the  {cq},  and  it  is  necessary  to 
resort  to  the  formulation  of  an  infinite  set  of  linear  equations  from 
which  (cq)  can  be  found  by  truncation. 


For  instance  in  the  case  a  -  0,  the  {c  }  must  be  such  that  the 

n 

equations 

f(1)  -  y  c  *(1)  (3-2a) 

“  n  n 

f(2)  -  l  c  *<2>  (3.2b) 

n  n 

are  simultaneously  satisfied  pointwiseon  0  g  y  i  1.  This  will  only  be 
possible  for  sufficiently  smooth  data  (f^l\  f^),  and  in  particular 


-9- 


if  Che  conditions 


f(1)(0)  -  0,  f(1)(l) 


0,  f  f  <2)(y)dy  -  0 


(3.3) 


(1)  (2} 

ate  satisfied  (since  each  of  the  $  7  f$  satisfy  the  corresponding 

n  n 

conditions).  Comparison  with  the  known  results  for  the  canonical  problems 

suggests  that  a  sufficient  condition,  for  data  satisfying  (3.3  ) ,  is  that 

(f  ^ ,  f^  belongs  to  the  Sobolev  space  H^(0,1),  i.e.  that  (f^  ,  f^  ) 
e  1*2  (0,1).  However  a  meaning  can  be  assigned  to  the  distribution  of  c^ 

obtainable  from  (3.2  )  under  considerably  more  general  conditions. 


Formulation  of  infinite  sets  of  equations  for  the  {cr  } 


Three  essentially  different  methods  for  formulating  equations  occur 
in  the  literature: 

(A)  Direct  collocation:  In  which  a  truncated  set  of  the  equations  0*2  a, b) 
are  identically  satisfied  at  a  suitable  set  of  points  {y  }  of  the  interval 

ID 


(o.i). 

(B)  Galerkin  methods  Here,  sets  of  weighting  functions  are 

—  ■  ■  —  ■  -  m  m 

introduced  and  mean  equations 

(W(D  f(l)}  .  £(w(D  *U>)  (W<2)  f(2))  -  I(W<2)^2))cn 

ID  IDQItDI  B  u  u 

are  obtained  and  solved  simultaneously.  The  simplest  such  equations  are 

provided  by  the  set 

■  sin  tmrx,  W^2^  ■  cosnnrx. 

m  m 


(3.4a,b, 


These  were  proposed  by  Benthem  (1963),  and  are  equivalent  to  those  used  by 
Gaydon  A  Shepherd.  With  our  notation  the  matrix  elements  are  given  by 


^in  nnry  ”  ^n^  .<^cos  mity  (~l)m  1 


2nnrX  tanX 
_ n _ n 

2  2.2 

B  IT  ) 


(3.5) 


( K 


(C)  Use  of  biorthogonal  functions 

Johnson  &  Little  in  effect  add  equations  (2.11)  and  (2.12)  to  obtain 


2c«  •<  1  >  <»:3,£<u  *  *w'<2,>  *  <|)  (^2>-  ^>fw> 

On  the  right  hand  side,  f^  and  f^  are  known,  while  f^  and  f^ 

are  calculated  as  £  c  T  c  <^4\  From  known  quadratures  we  find 

n  n  *•  n  n 


(3.6) 


the  resulting  bracket  as 

1  ('♦«)♦<»♦(♦»).  «<4>)*<4A 

2  \Hm  yn  vym  *ra  /Hn  / 

say,  while  its  value  is  1  for  m  *  n. 

-11- 


■2^_(^_tanX  -X  tanX  ) 
mm  m  n  n 


F  (mj»n) 
mn 


(X  -X  )^(X  +X  ) 
m  n  m  n 


>(3.7) 


i 


This  gives  the  set  of  equations 


c  »  7  P  c  ♦  d 
n  1  n  n  a 


where 


(3.8) 


-  \  ♦r)‘<2)> 


It  is  to  be  noted  however  that  this  set  of  equations  can  also  be 
obtained  directly  by  a  Galerkin  method  from  Q. 4a.  and  b),  using  the 
weights  V<1}  -  \  «£3).  -  \  adding . 


For  an  infinite  set  of  equations  of  the  form  (  3.8),  Kantorovich  &  Krylov  (1958) 
show  that  the  solution  {c  }  is  stable  in  the  sense  of  approaching  a  limit  as 

D 

the  number  N  of  equations  become  large  if  the  L  norm  li?||  =  sup  I  |P  I  is  <  1~£,£>0. 

oo  oo  m^n  *  mn 

However  this  condition  is  not  satisfied  by  the  coefficients  (3.7  );  in  fact  the 

sum  y  |F  |  tends  to  infinity  with  m  . 
n  0,1 

To  show  this  the  equations  may  he  written 


c 

m 


(F, 


mn 


c 

n 


+  F 


mn 


Cn> 


+  d 


where 


mn 


-  2A  (A  tanA  -X  tanA  )/(A  -X  )2(X  +A) 
n  u  mn  n  m  n  mn 


(3.9) 


The  summation  now  extends  over  Xr  in  the  first  quadrant  only,  and  we  need 

only  consider  X  in  the  first  quadrant  since  the  conjugate  equation  holds  for 

m 

A  .  The  asymptotic  expressions  (l.ll)  show  that 

D 

X  ,  X  nm,  while  tanX  ,  tan  X  -*•  i. 
m  o  m  n 


Therefore 


m 


(m-n) 


-2 


(m  +  n) 


(3.10) 


-12- 


*  2  * 

(For  m  ■  n,  ■  0,  F  A#  2inm/(£n4mir)  ).  The  imaginary  part  of  F  does 
on  mn  mo 

not  appear  in  the  diagonal  when  the  equations  (3.8  )  are  formulated  in  real 
and  imaginary  form.  Then 


^  ~  in  2m,  £  |F*  J 


nnr 

3 


(3. 


Similar  behaviour,  ia  exhibited  by  the  elements  (3.5  )  of  the  matrix 
formed  by  sine  and  cosine  weighting  functions,  and  we  might  expect  this  to 
be  reflected  in  instability  of  the  solution  of  truncated  sets  of  these 
equations.  In  fact  this  was  shown  in  numerical  experiments  described  infl]. 


It  is  evident  that  the  difficulties  arise  from  the  factor  (A  -A  )*  in 

m  n 

the  denominator  of  F 

on 

Some  light  Is  thrown  on  the  behaviour  of  the  Inverse  matrix  (I  -  F) 
by  the  distribution  of  eigenvalues  of  I  -  F.  For  N  ■  4,  l.e.  the  8^8 
matrix,  the  eigenvalues  are 

I  -  F  (Johnson-Little):  1  ±  1.3749,  1  t  .8337,  1  ±  .9901,  1  ±  1.0005  . 

Thus  the  last  2  eigenvalues  are  extremely  close  to  zero,  and  the  matrix 
Ill-conditioned  (det  ■  5  *  10  ®) . 

By  contrast,  for  the  matrix  I  -  G  derived  In  the  next  section  by  use 
of  optimal  weighting  functions,  the  eigenvalues  are 

I  -  G:  1  ±  .3920,  1  i  .0828,  1  ±  .0060,  1  1  .0001 


and  It  appears  that  the  eigenvalues  rapidly  approach  1  for  larger  matrices, 
with  det(I  -  G)  of  order  1. 


Optimal  weighting  functions 


As  noted  in  the  last  section,  the  matrices  resulting  from  the 
Galerkin  methods  B  and  C  are  not  diagonally  dominated,  and  this  fact 
is  reflected  in  the  computations  described  in[  1  ],  where  it 
was  found  that  the  coefficients  obtained  from  these  methods  varied  with 
the  order  of  truncation. 

However,  an  extra  degree  of  freedom  not  exploited  by  previous 
authors  is  available  in  the  choice  of  weighting  functions.  In  the  hope 
of  producing  a  diagonally-dominated  matrix,  consider  weighting  functions 
of  the  general  form 

-  * ♦«'  . •  c 

-  c  ♦<»  .  0  *<‘> 


This  is  a  natural  extension  of  the  biorthogonal  weighting  functions  (3.7  ), 

and  contains  sufficient  disposable  constants  to  ensure  the  suppression  of  the 
2 

factor  (X  -  X  )  in  the  denominator  of  the  matrix  elements, 
m  n 

If  we  write 


p(D 


♦  a  f 


O)  .  g(U. 


e<2> 


-  a 


f(4>  .  g(2) 


(3.1  bis) 


then  the  stress  (a  -  0)  and  displacement  (a  ■  1  -  v)  problems 

together  as  cases  in  which  the  functions  g^  are  known 

shall  therefore  seek  optimal  weighting  functions  in  the  above 
(2) 

and  g  ,  and  expect  the  result  to  depend  on  the  parameter  a. 
equations 


.(1) 


z  =„  <♦»>  ♦  «♦»>> 


.«> 


«  '  ■  Z  c„  <*‘2)  -  a*™) 


.(4), 


may  be  treated 

on  CKy<l .  We 

sense  for  g^ 
Thus  from  the 


(3.13) 
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ditt 


we  obtain  the  equations 


y  A _  c  •  d 

*•  an  n  n 


a„  .  *"vw— :*)> 

A.  Xi2>.<2>> 

Quadratures  of  the  eigenfunctions  appendix  of  (1)  show  that  for  m  ^  n 


(X2-X2) 
n  n 


A(>i  •  A2)  « 

X2  “  X2  \  XX 

An  n  \  m  n 


where  C 


♦  4  a|(B-D)  C^-BT^I 

m/c°B\my  cosXny\  .  /x  tan  X  -  3 

\  cosX  cosX  /  \  m 

\  n  n  ' 


X  tan  X  -  A  tan  A  )/(A2  -  X2  ) 

n  tn  n  n  m  u 


and  T 


(tan  A_  ♦  tan  XJ  /  (X^  ♦  Xn) 


(3.14) 


(3.15) 


while 


A  •  iA  ♦  |  C  ♦  |  (X  tan  X  )  (A  ♦  C)  ♦  B  ♦  D 
ma  2  3  a  n 


♦  a  I  A-C  -  4 


■&■)] 


To  suppress  the  factor  (Xn  -  Xn>3  in  the  denominator  of  it  is  necessary  to 
set  A  -  -  C.  The  choice  C  -  1,  A  -  -1  reduces  the  expression  to 


'  XX 

a  n 

<w: 


(B+o)X  X  ♦  X2  D  ♦  X2o  | 

- 5LD - 2 - S-  ♦  a(B-D) 

X2  -  X2  J 


C  -4aBT__ 
mi  uui 


a  n 
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If  now  wo  sot  B-J-a,  D--J-a,  chc  remaining  factors  of 
in  the  denominator  are  suppressed  leaving 


A _  »  -  G _  -  2a  (i-2d)  T  m  +  n 

Bu  in  on 


•ay.  where  G _  -  2X  (X  tanX  -  X  tanX  )  /  (X  +X  )' 

on  o  o  on  n  on 


while 


(l-4a)  -  2a(l-2a)T 


> 


Hence  if  the  weighting  functions  are  chosen  as 


♦  <l-ot)  ♦„3) 


♦iZ)  “  <i+«>  *i4) 


the  equations  (3.14)  are  of  the  form 


(l-4a)  c  -  l  (  G_  ♦  2a  (l-2a)T  >  c  ♦  d 

a  on  an  n  a 


In  particular,  for  the  traction  problem  (a  ■  0)  they  reduce  to 


c  ■  F  G _  c  ♦  d 

m  **  mn  n  m 


Likewise  for  the  Stokes  problem  (a  »  1)  the  system  is 


-c  ■  y  G  c  ♦  d 
m  u  am  n  a 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


(3.20) 
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Diagonal  dominance  of  the  matrix  A  5  1  -  G 


We  now  confirm  that  the  infinite  matrix  I-G  is  strictly  diagonally 
dominated,  i.e.  that 

l  |0ml  <  1  ».  <3-21> 

n— « 

An  estimate  of  the  absolute  row  sums  of  G  may  be  obtained  by  use  of 

mn 

the  asymptotic  values 

».  *  *«  '  <l«l  -  £> « 

X.  t-X.-ijli-  X*)‘  <.-1;  IX.  (.  ?  0) 

Hence  in  the  numerator  of  I G _ I , 

mn 

|AmtanAm  -  X^tanXj  (m  -  ^  )  ♦  (n  -  j  )  |ir 

accordingly  as  m,  n  have  the  same  or  opposite  signs. 

Likewise  the  denominator  of  |G  I  namely 

mn 

|Xb  +  Xn|3  >  (ReXa  +  ReXn)3  %  it2  (|m|  +  |n|-  |)3 
Hence  for  m,  n  >  0, 

'(*)  [  l<Siqft!->)]  5  «<■•">  <*•*» 

This  leads  to  an  estimate  of  the  mC^  absolute  row  sum  of  the  truncated  matrix 
of  order  2N  as 

H  N  /.N+i 

l  |6j  *  l  g(»,n)  <  g(m,n)  dn  (3.23) 

n«-N  n-1  J  | 

(the  inequality  holds  since  g  is  decreasing  and  concave  as  n  increases). 
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In  particular,  in  either  case  d  •  0  c  *  0 


(3.30) 


Table  2 


m 

30 

l  |CJ 

n— 30 

S  -R  (3°> 

■  m 

99 

J  i°J 

n— 99 

s  -  »  <”> 

n  b 

1 

.63515 

.64598 

.65744 

.66690 

2 

.74430 

.73800 

.79091 

.78385 

3 

.75863 

.75157 

.82716 

.81938 

4 

.75142 

.74515 

.83923 

.83233 

5 

.73713 

.73179 

.84194 

.83603 

6 

.72040 

.71586 

.84021 

.83517 

7 

.70306 

.69918 

.83610 

.83177 

8 

.68593 

.68257 

.83063 

.82688 

9 

.66937 

.66644 

.82436 

.82108 

10 

.65355 

.65097 

.81759 

.81470 

15 

.58494 

.58646 

.78152 

.77980 

20 

.53652 

.53551 

.74648 

.74534 

25 

.49894 

.49822 

.71432 

.71349 

30 

.47011 

.46956 

.68522 

.68459 

40 

.63535 

.63495 

50 

.59471 

.59443 

60 

.56128 

.56107 

70 

.53350 

.53334 

80 

.51019 

.51006 

90 

.49044 

.49034 

99 

.47515 

.47506 

The  first  and  second  columns  are  the  two  sides  of  (3.26)  for  N  -  30, 
the  third  and  fourth  for  N  ■  99. 
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4.  Asymptotic  estimates  of  convergence  retes 


In  this  section  asymptotic  estimates  vill  be  made  of  the  rate  of 
convergence  of  the  coefficients  {cq}  for  the  mixed  problem  in  which 
f^,  f^  are  prescribed  functions  of  y  on  the  interval  [0,1]. 

Four  types  of  data  will  be  considered 

(i)  Smooth  continuous  data  satisfying  the  comparability 
condition  f^(l)  -  0. 

(ii)  Data  for  which  the  above  comparability  condition  is 
violated. 

(iii)  Data  containing  a  discontinuity  at  an  internal  point 
c  e  (0,1). 

(iv)  Data  containing  a  discontinuity  in  slope  at  an  internal  point  c. 

In  each  case  it  will  be  assumed  that  the  data  satisfies  the  additional 
requirements  at  y  ■  0  for  expansion  in  terms  of  the  even  eigenfunctions. 

These  are 

f(1)(0)  -  f(1)"(0)  -  0,  f<2)'(0)  -  f<2)"(0)  -  0.  (*•!) 

The  technique  in  cases  (i)  and  (ii)  will  be  to  estimate  the  behaviour  of 
the  functions  f^and  which  is  not  known  in  advance,  by  solving  the 

bihaxmonic  equation  in  a  quarter  plane  with  boundary  conditions  on  the  two 
sides  that  are  asymptotically  the  same  as  those  in  the  neighbourhood  of  the 
corner  x  ■  0,  y  »  -  l  in  the  strip  problem.  In  cases  fcii)  and  (Lv),in  which 
the  governing  singularity  in  the  data  is  at  an  interior  point,  the  relevant 
asymptotic  solution  of  the  biharmonic  equation  is  in  a  half-plane.  Having 
obtained  estimates  of  f^  and  f^\  we  can  then  estimate  c  from  (2.11)  or 
(2.12),  using  in  addition  the  known  values  of  f^  or  f^2\ 


where  Y  »  1  ♦  y. 

The  solution  of  the  bihamonic  equation  in  the  quarter  plane  X,  Y  £  0, 


The  value*  of  the  Laplacian  P  -  AT  and  its  harmonic  conjugate  Q  on  the 
boundary  6  ■  */2  are  obtained  as  the  inverse  Mellin  transforms 


where 


(:)  •*!'(■:)■•“ 


-  ■(■-) 


z  ^  .  2  ITS  .  .  ITS  _  _  ITS 

(i+nn  y  )  sxn  y  cos  y 


1f|  ITS  g  % 

•IQ  Y~  cos  2^  (•  “  •in*^~) 


and  A  ■  det  M  -  a*  -  sin2  y 


(4.6) 


(4.7) 


A  has  real  seros  at  0,  0,  ♦  1,  and  a  double-infinite  set  of  complex  zeros  ♦  y2, 
+  y  ,  ♦  Y.»  *  Y a  » • • •  which  are  spaced  so  that  2n  <  ReY2n  <  2n  ♦  1 


The  first  of  these  is 


Y,  -  2.739593  ♦  i  (1.119025) 


(4.8) 


The  Mellin  transforms  of  g  and  h  must  be  of  the  form 
^  *2 

h  (s)  ■  -  — j  “  "J73  +  function  regular  in  -  3  <,  Res  <0  J 

*  k  b.  b.  J 

g(s)  *  — ^  ♦  — j  ♦  — j  ♦  function  regular  in  -3  <  Res<  q 

since  on  inversion  in  accordance  with  (  4.5),  these  produce  the  stated 
expansions.  Then  from  (4.7), 


■  (-  !)(•♦•*”'¥)  *  *  (1xa) 


p  (*) 


(4.10) 


On  inversion  using  (4.6)  we  obtain 


p(Y)  -  b  ♦  b.Y  ♦  CYY2“l  +  C  YYj“1  ♦  0(Ya) 
o  i 

where  C  is  a  (complex)  constant. 


(4.11) 


Since  f 
interval 


<2).  f 

0  <  y 

f<2> 

f(4> 


<4)  are  sytanetrical  with  respect  to  y,  they  can  now  be  expressed  on  the 
<  1  by  writing  1-y  for  Y  in  (4.2b)  and  (4.11)  respectively,  as 

-  bo  ♦  bj  (1-y)  +  0  (1-y)2 

-  b  ♦  b.  (1-y)  ♦  C  (l-y)Ya  “  1  ♦  C  (l-y)Yt  “  1  ♦  ... 

o  i 


(4.12a) 

(4.12b) 


c  can  now  be  evaluated  from  (2.12)  as 


c 

n 


i  <-  (*s£)  «<2>  -  ‘w> 


Two  integrations  by  parts  give 

1  /  2cosAny 

n  2A  2  \  cosA 

n  \  n 


2cosA  y 
n'' 

A2cosA 
n  n 


(4.13) 


For  f^2\  f^  C  H2  (0,1),  bounds  on  these  terms  are  provided  by  the 
Cauchy-Schwarz  inequality,  since  from  [1]  we  have  the  norms 

*>»  const  /  (In  n) ^ 


2cosAny 

cbsX_ 


(4.14a 


||  ||  ^  const  /  n  (In  n)2<^2 

At  this  stage  therefore  our  best  estimate  is 


(4.14b‘ 


I Cn  |  •  0  (l/n  (In  n )  ^  , 

(4) 

since  f  cannot  be  expected  to  have  stronger  differentiability  properties 
than  the  primary  data. 


(4.15) 


i 


(2)  4 

However  for  smoother  data  such  that  f  e  H  (0,1),  we  may  write 


(2)" 

-  u(y) 

► 

(4.16) 

(4)" 

-  K  (l-y)V_1  ♦  K  (l-y)V”1  +  v(y) 

J 

where  K  -  (y  -1)(Y  -2)C,  v  ■  y  -2,  and  u(y),  v(y)  are  such  that  u",  v"  e  L_  (0,1) 

2  2  2  t 


Then 


'4 


2cosX  y 
_ n 

X2cosX 
n  n 


(u  -  v)  +  ifiQ  v 


<(• 


2cosX  y 

+  i<l*n  _  ~~i — —  II K  (1_y)V  A  +  K  (1 


X  cosX 
n  n  ■ 


-•J> 


(4.17) 


”3. 


Further  integration  by  parts  shows  that  the  first  term  is  0(n  ) . 

We  now  focus  attention  on  the  second  bracket.  The  part  multiplied  by 
K  when  written  in  full  is 


(2Acos2A) 


1  /  r  • 

(  {si 


inX(l-y)  ♦  (1-y)  cosXsinXy  -  t-  cosXcosXy}(l 


X-X 


(4.18) 


After  integration  of  the  middle  term  by  parts,  and  on  writing  y  -  1-t,  this 
becomes 


(2Acos2X)-1  ^{sinXt  -  (2+v)  cosX  (1-t)  }  t 


+  (2X2cosX)_1 


(4.19) 


The  integral  can  be  estimated  for  large  n,  by  use  of  asymptotic  formulae  derived 
in  Appendix  D  for  X  e  A  ♦,  i.e.  the  set  of  X  for  which  C(A)»  0,  ReX  >  0,  and 
0  <  Rev  <  1: 


<■ 


tV  *  cosX^  ■  X  v(cos  —  )  f(v)  ♦ 


<^V-X  sinXt^  ‘  ^ 


/ 


X  (sin  )  r(v)  - 


sinX 

X 

-3/2 

♦  0(A  ) 

(4.20a) 

cosX 

X 

-  3/2 

+  0(X  ) 

(4.20b) 
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When  theaa  ir«  aubatituted  in  (4.19)  th«  contribution*  iron  tho  aecond 

2 

Coma  on  tho  right  cancel  perfectly  with  the  tern  1/2X  co*A,  end  we 
are  left  with 

(’  l)^"2  r<v>  £<**v>  y  ♦  O  ♦  v)  (Bin  y  )tanxj  ♦  0(X  ?/2> 

Since  v  •  Y«“  2*  where  Y.ia  auch  that  Y  ♦  ain  0, 

2  2  2  « 

thia  expreaaion  can  be  written  in  the  form 

x”*2(a  4  |  tanX)  ♦  0  (X  ^2)  1 


(4.21) 


Then  aince  X  %  nw  ♦  0  (■“—).  while  tanX  ^  ♦  i,  the  estivate  of 


leBl  i* 


|c  |  "viconaOn  ”2,7*  ♦  0(n"3) 


Thia  estivate  ia  well  aupported  by  coaputatione  using  the  aethod  of 
optiael  weighting  function*,  for  two  caaea  for  which  the  coefficient*  cR 

were  tabulated  in  [1].  In  the  noaenclature  of  that  report,  they  are 


Ceae  1«  f(l)  •  0,  f<2> 


Ceae  )i  f 


y  -  y*.  f<2> 


(4.22) 


(4.23) 


The  value*  of  n  '  *  ^  ,  with  H  •  40,  are  liated  in  Table  2  ,  and  it  ia 

aeon  that  in  both  caaea,  they  vary  only  alowly  with  n,  the  vaxivuv  variation 
between  n  •  12  end  n  ■  40  being  of  order  ♦  8  per  cent. 


Tabic  3  Decay  of  coefficients  for  the  smooth  continuous  data  listed  in  (4.23) 


Case  1 

Case  3 

Case  1 

Case  3 

.  n 

2.74,  , 

n  lcJ 

2.74,  , 

°  lCJ 

n 

2.74.  , 

n  lcnl 

2.74,  , 

»  |cnl 

1 

.8676 

.3704 

21 

.4555 

.4854 

2 

.4735 

.6337 

22 

.4566 

.4945 

3 

.6175 

.6465 

23 

.4586 

.5030 

4 

.6774 

.6012 

24 

.4615 

.5108 

5 

.6912 

.5473 

25 

.4656 

.5179 

6 

.6820 

.4978 

26 

.4697 

.5243 

7 

.6616 

.4586 

27 

.4742 

.5301 

8 

.6362 

.4291 

28 

.4790 

.5352 

9 

.6093 

.4096 

9 

.4839 

.5398 

10 

.5830 

.3988 

30 

.4886 

,  .5438 

11 

.5583 

.3951 

31 

.4940 

.5473 

12 

.5361 

.3969 

32 

.4986 

.5503 

13 

.5167 

.4027 

33 

.5034 

.5529 

14 

•  50c 

.4111 

34 

.5078 

.5550 

15 

.4865 

.4211 

35 

.5120 

.5567 

16 

.4756 

.4320 

36 

.5164 

.5580 

17 

.4675 

.4432 

37 

.5202 

.5590 

18 

.4615 

.4544 

38 

.5233 

.5597 

19 

.4578 

.4652 

39 

.5270 

.5600 

20 

.4560 

.4756 

40 

.5295 

.5601 

Note: 

The  calculation  for  Ca&e  3 

was  also  repeated  with  N  *  99, 

and  values  of 

Y, 

|X  2c 

1  showed 

similar  behaviour 

>  e.g. 

n' 

n 

20  30 

40  50 

60  70 

80  90  99 

Yo 

lXn  c 

|xlO~2 

n‘ 

.10192  .12137 

.12899  .13002 

.12777  .12407 

.11998  .11611  .11308 

with  smooth  variation  between  these  values.  The  values  differ  from  those  in  the  table 

y 

by  the  lector  |An  2|/n2’^  which  ^  ir2,2*  -  23.0245  as  n  -*■  ®  . 
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(ii)  Data  not  satisfying  the  edge  condition  f^(l)  “  0 


To  investigate  the  effect  of  data  that  violate  the  compatability 

condition  4*  (0,1>  -  0,  it  is  sufficient  to  consider  the  case 

xy 

f(l)(y)  -  y.  f<2)(y)«0 

since  smooth  data  of  the  type  already  considered  can  be  superposed  on 
this  to  produce  a  general  distribution.  Referred  to  the  corner  (0,-1), 
f(1)  •  -  1  +  Y  so  that 


(4.24) 


;a> .  .  1  +  1  + 

S+l  8+2 


~  *(1)  '■ 

In  this  case  from  (4.7)  with  h  ■  f  ,  g  ■  0,  we  have 


Csinirs  \ 
4  ) 


<i) 


which  has  a  simple  pole  at  s  «  -1  with  residue  y  and  is  regular  at  s  •  -2, 


the  next  poles  being  those  at  -y  ,  -y 

2  2 


Therefore 


.-I* 

and  leading  term  in  cr  is  obtained  by  writing  f^  *  y,  f^  «  0  in  (2. 12). This  gives 


e  *  -y-  .  y/  -X2  \ 

n  2  2\  n  “  ^ 


TT 

2 


rtanX 


(4.25) 


—  Y 

the  next  terms  being  0(X  2)* 

n 


-27- 


■  i  - 


Since  can  X 

II 


-2 

♦  0  (X  ),  this  result  shows  Chat 


|c  X  I  ■»  ?  «s  n  ♦  •  (4. 

n  n  i 

which  is  well  borne  out  by  a  calculation  using  the  method  of  optimal 

weighting  functions,  as  shown  in  table  4.  Except  for  n  ■  1,  the  computed 

value  of  IX  c  I  lies  within  2  percent  of  y  -  1.5708,  and  the  coefficients 
1  n  n 1  ^ 

in  fact  approach  the  pure  imaginary  fcrm  Jn 


implied  by  (4.25) 


Table  4  Data  violating  tha  compatibility  condition 

3tanl 

fH)  .  yf  f(2)  •  o  d^  -  -  y  +  - 2 — 

X. 

The  coefficient*  were  computed  with  H  -  99. 


b 

|A  c  | 

n 

*n 

“n 

n  n 

1 

0.243986, +00 

0.777206, +00 

0.194527 ,+01 

2 

0.376416,-01 

0.277808, +00 

0.1 56334, +01 

3 

0.213641 ,-01 

0.174970, +00 

0.1 53696 ,+01 

4 

0.142829,-01 

0.1 29271, +00 

0.154212, +01 

5 

0.101696,-01 

0.102807 ,+00 

0.154905,+01 

6 

0.755116,-02 

0.853834,-01 

0.155421  ,+01 

7 

0.578495,-02 

0.730035,-01 

0.155763,+01 

8 

0.453904,-02 

0.637432,-01 

0.155979, +01 

9 

0.362747,-02 

0.565529,-01 

0.1561 05, +01 

10 

0.294002,-02 

0.508079,-01 

0.156168, +01 

11 

0.240808,-02 

0.461126,-01 

0.1 56188 ,+01 

12 

0.198738,-02 

0.422037,-01 

0.156176, +01 

13 

0.164834,-02 

0.388992,-01 

0. 156142,+01 

14 

0.137055,-02 

0.360694,-01 

0.1 56091 ,+01 

15 

0.113968,-02 

0.336189,-01 

0.156027 ,+01 

16 

0.945324,-03 

0.314764,-01 

0.1 55954, +01 

17 

0.779874,-03 

0.295875,-01 

0.155874.+01 

18 

0.637593,-03 

0.279097,-01 

0.155789, +01 

19 

0.514141,-03 

0.264095,-01 

0.1 55700, +01 

20 

0.406145,-03 

0.250602,-01 

0.155608, +01 

21 

0.310976,-03 

0.238402,-01 

0.155514,+01 

22 

0.226564,-03 

0.227318,-01 

0.155418, +01 

23 

0.151235,-03 

0.217204,-01 

0.155320 ,+01 

24 

0.836376,-04 

0.207937,-01 

0.1 55223 ,+01 

25 

0.226839,-04 

0.199417,-01 

0.1 551 24 ,+01 

26 

-0.325382,-04 

0.191555,-01 

0.155025, +01 

27 

-0.827604,-04 

0.184280,-01 

0.154927 ,+01 

28 

-0.128630,-03 

0.177527,-01 

0.1548 28, +01 

29 

-0.170660,-03 

0.171243,-01 

0.154730 ,+01 

30 

-0.209291,-03 

0.165381,-01 

0.154633, +01 

31 

-0.244919,-03 

0.159899,-01 

0.1545 36, +01 

32 

-0.277862,-03 

0.154762,-01 

0.154439, +01 

33 

-0.308389,-03 

0.149938,-01 

0. 154344, +01 

34 

-0.336752,-03 

0.145400,-01 

0.1 54249, +01 

35 

-0.363157,-03 

0.141122,-01 

0.1S4155,+01 

36 

-0.387793,-03 

0.137084,-01 

0.154062, +01 

37 

-0.410818,-03 

0.133266,-01 

0.153970, +01 

38 

-0.432371,-03 

0.129649,-01 

0.153878, +01 

39 

-0.452582,-03 

0.126220,-01 

0.153788, +01 

40 

-0.471558,-03 

0.122963,-01 

0.153699, +01 
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(iii) 


Data  containing  discontinuities 


We  next  consider  the  case  when  one  of  the  date  functions  f^,  f*'*' 
has  a  finite  discontinuity  et  an  internal  point  y  ■  c  of  (0,1),  and  is 
elsewhere  continuously  differentiable.  In  this  case  the  solution 
sufficiently  close  to  the  point  c  is  in  the  limit  the  same  as  in  a 
half  space  x  >  0.  For  the  half  space,  the  following  singular  integral 
relation  exists  between  the  second  derivatives  at  the  surface  x  ■  0  of  a 
function  V  satisfying  the  biharmonic  equation  (appendix  E  ) 

-  v  <«•*  -  ?  £  y°-cuc . 

the  integral  on  the  right  being  a  Cauchy  principal  value.  The  integral 
is  continuous  at  a  point  of  continuity  of  ¥^(0,t),  and  is  continuously 

differentiable  if  ¥  (0,t)  is  so.  Therefore  a  point  of  continuity  of  the 

xy 

boundary  value  of  ¥  is  also  a  point  of  continuity  of  ¥  -  ¥  ,  and  if 

xy  xx  yy 

¥^  has  a  jump  discontinuity  at  a  point  y  -  c  of  the  boundary,  has  an 
equal  jump  at  the  same  point,  and  A¥  a  jump  of  twice  the  amount. 


Thus  for  the  strip,  if  f^  is  continuous  at  c  while  f^  has  a 
discontinuity  of  amount 


f(2\c+)  -  f(2)(c-)  =  [f(2)1  say 


(n 

there  will  also  be  a  discontinuity  in  f  of  amount 


f^(c+)  -  f(4)(c-)  .  2  [f(2^] 


(4.27) 


(4.28a) 


(4.28b) 
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Th«  contribution  to  c  is  therefore  given  by  (2.12)  es 

n 

iff  f<2)<c->  ♦n4>  ♦  <(*W2>-  **4)>  J  dy 


iff 


t(2U  ♦«>.  f 


‘W’-  C) 


(2)  " 

The  dominant  contribution  is  that  from  $  =  to  ,  and  is 

Q  u 


t-i)  u‘*w  *»vi  •:» 


/  *ln>ne  .  “•  V  \ 

1  5c  \  ,io>n  ‘  ) 


■v  -  tf<2)]c  <1-0  «ixa(l”c) 


Since  Im  ^  j  In  4nir,  we  obtain  from  this  expression  the  estimate 

•  (il£) 

|cn|  ^  (1-c)  |{f(2)]c|(4nir)  2 


In  [1] ,  the  coefficients  generated  by  the  distribution 

,(2> 


f(l)-  0. 


f  1  0<y<  | 

l-i  i  <  y  <  i 


were  calculated  using  the  method  of  optimal  weighting  functions.  This 

falls  into  the  pattern  just  discussed,  with  c  -  j,  [f^l  *  -  2.  The 

c 

estimate  for  this  case  is  therefore 

1/4  ,  ,  -1/4 

n  | cn |  "*  (4w)  •  0.5311 


(4.29) 

dy 


(4.30) 

(4.31) 


(4.32) 


This  is  borne  out  very  well  by  the  calculated  values,  as  shown  in  Table  S  . 
Alternate  values  are  above  and  below  .5311,  the  means  of  successive  pairs  of 
values  showing  rapid  convergence  towards  this  value. 


Case  2:  « 


0, 


:<2) 


1  y  £  (0,J*) 
-lye  (%,1) 


For  this  case  d 


tan 


X 

_  sin 

-  +  2+  '  Vi)  +  2  - - Lil  ‘ 


m 

(The  last  two  terns  can  be  combined  as 


X  cosX 


tanX 

m 

i  f3X  ) 

1  ml 

,  3 

M 

X 

m 

-  j  c°8(-2-J 

+  2 

cos 

[Tj 

the  form  given  In  [1]). 

C»  "  +  lbn  '“P"'"1  *lth  »  -  40  nr.  IK  ted  on  page 


40 


1/4.  , 

"  lcJ 

n 

1/4,  . 

"  u 

Mean  of  successive 

.86252 

21 

.60434 

values 

.49262 

22 

.47357 

.5389 

.73997 

23 

.60258 

.5380 

.46516 

24 

.47255 

.5376 

.70536 

• 

25 

.59771 

.5351 

.43944 

26 

.47743 

.5376 

.56952 

27 

.59695 

.5372 

•45011 

28 

.47812 

.5365 

.65088 

29 

.59260 

.5354 

.45364 

30 

.48059 

.5366 

.63447 

31 

.59243 

.5365 

.45633 

32 

.47912 

.5358 

.62715 

33 

.58852 

.5338 

.46245 

34 

.48325 

•  5359 

.61982 

35 

.58870 

.5360 

.46286 

36 

.48168 

.5352 

.61347 

37 

.58516 

.5334 

.46875 

38 

.48552 

.5353 

.60988 

39 

.58556 

.5355 

.46821 

40 

.48389 

.5347 

(tv)1'*  - 


0.531126 


(iv) 


Discontinuity  in  a  derivative 


Continuing  the  study  of  discontinuous  data,  suppose  next  that  f 
is  continuous,  but  has  s  jump  in  first  derivative,  of  amount 


(2) 


f(2)(c+)  -  f(2)(c-)  -  [f(2)'l 


(4.33a) 


at  a  point  c  e  (0,1).  Then  from  the  last  section  it  follows  that  provided 
f^  is  continuously  differentiable  at  c,  there  is  a  corresponding  jump  in 
f  ^  ,  of  amount 

f(4)(c+)-  f(4)(c->-  2  [f(2)']c  (4.33b) 

The  expression  (2.12)  £or  cq,  after  one  integration  by  parts,  gives 

%  -(-?)< 

n 

(where  we  have  also  used  the  fact  that,  for  compatible  data, 

f(4U)  -  f<2h)) 

The  dominant  term  is  the  last  in  the  bracket,  which  contributes 

cn  (f<4)(c+)-  f<4)(cr^n(c)  (4*34) 

As  before,  asymptotic  results  show  that 

V*> '  (in) 

whence 

I'J  '  (sr) 


iA  (i-c)  . 

(l-c)e  ,  X  *v»  nir  +  •=■  In  Ann 

n  i. 


(407T)  (  2  }  |[f<2)’]c| 


(4.35) 
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The  coefficients  generated  by  the  distribution 


fd) 


•  0, 


r  1  0  <  y  <  J 

1  -  8(y-i)  4  <  y  <  i 

> 


(4.36) 


have  also  been  calculated  by  the  method  of  optimal  weighting 

functions,  and  are  listed  in  Table  6.  For  this  case  c  •  4, 
/ 2 1 1 

[f  1  ■  -  8  so  our  estimate  is 

c 


(nir) 


5/4 


-5/4 

.67625n 


(4.37) 


and  this  is  borne  out  by  the  computations,  again  with  oscillations  above  and 
below  in  successive  terms. 


The  right  hand  side,  calculated  from  (4.36),  is 


ij;"  + 


3tanA_ 


-  0) 


8 

F 

m 


(4.38) 


The  second  term  here  equals  -  8<j/(|),  and  has  precisely  the  same  asymptotic 
behaviour  as  (4.37). 


T»bl>  6  Computations  made  using  the  method  of  optimal  weighting  functions 
RHS  is  (0,  1  0  <  Y  <  0.5}  5  -  8y  0.5  <  V  <  1) 

Number  of  equations  is  50 


Mean  of 


*n 

bn 

»5/4i=»i 

successive 

c 

1 

0.85217905, +00 

0.82889897 ,+00 

1.18889 

c 

2 

0.15110750, +00 

0.3127 3898, +00 

.82610 

c 

3 

-0.23383817 ,+00 

-0.92803374,-01 

.99330 

c 

4 

-0.1 1630692, +00 

-0.1 7379686, +00 

1.18298 

c 

5 

0.81507437,-01 

-0.16184736,-01 

.62131 

c 

6 

0.53842954,-01 

0.69301344,-01 

.82411 

c 

7 

-0.52680854,-01 

-0.26550267,-02 

.60059 

c 

8 

-0.33598252,-01 

-0.57684090,-01 

.89815 

c 

9 

0.39962947,-01 

-0.71828929,-02 

.63294 

c 

10 

0.26417607,-01 

0.38246136,-01 

.82660 

c 

11 

-0.28607324,-01 

0.36837519,-02 

.57782 

c 

12 

-0.18216213,-01 

-0.31931769,-01 

.82107 

c 

13 

0.24928523,-01 

-0.48384923,-02 

.62684 

c 

14 

0.16112214,-01 

0.25319424,-01 

.81273 

c 

15 

-0.19434066,-01 

0.40120482,-02 

.58579 

c 

16 

-0.12109066,-01 

-0.21474948,-01 

.78892 

c 

17 

0.17645042,-01 

-0.37712050,-02 

.62285 

c 

18 

0.11117389,-01 

0.18456077,-01 

.79883 

c 

19 

-0.14567490,-01 

0.35850678,-02 

.59511 

c 

20 

-0.88840248,-02 

-0.15939817,-01 

.77181 

c 

21 

0.13460482,-01 

-0. 312131 S5, -02 

.62117 

.69649 

c 

22 

0.82831224,-02 

0.14294083,-01 

.78714 

.70416 

c 

23 

-0.11554256,-01 

0.31138903,-02 

.60274 

.69494 

c 

24 

-0.69142859,-02 

-0.12550230,-01 

.76116 

.68195 

c 

25 

0.10782529,-01 

-0.26671130,-02 

.62093 

.69105 

c 

26 

0.64976760 ,-02 

0.11540997 ,-01 

.77759 

.69926 

c 

27 

-0.95120892,-02 

0.27102167,-02 

.60874 

.69317 

c 

28 

-0.55983579,-02 

-0.10276338,-01 

.75374 

.68124 

c 

29 

0.89361062,-02 

-0.23263043,-02 

.62142 

.68758 

c 

30 

0.52875285,-02 

0.96042473,-02 

.76976 

.69559 

c 

31 

-0.80422256,-02 

0.23802368,-02 

.61350 

.69163 

c 

32 

-0.46653040,-02 

-0.86527597,-02 

.74818 

.68084 

c 

33 

0.75939475,-02 

-0.20592124,-02 

.62232 

.68525 

c 

34 

0.44224351,-02 

0.81778625,-02 

.76330 

.69281 

c 

35 

-0.69375644,-02 

0.21117078,-02 

.61735 

.69033 

c 

36 

-0.39725085,-02 

-0.74400338,-02 

.74374 

.68055 

c 

37 

0.65785572,-02 

-0.18437055,-02 

.62345 

.68360 

c 

38 

0.37777461,-02 

0.70893951 ,-02 

.75790 

.69068 

c 

39 

-0.60795518,-02 

0.18915114,-02 

.62053 

.68922 

c 

40 

-0.34414227,-02 

-0.65025212,-02 

.74008 

.68031 

c 

41 

0.57856459,-02 

-0.16661447,-02 

.62464 

.68236 

c 

42 

0.32814420,-02 

0.62349685,-02 

.75334 

.68899 

c 

43 

-0.53957403,-02 

0.17088882,-02 

.62322 

.68828 

c 

44 

-0.30225881,-02 

-0.57580215,-02 

.73695 

.68009 

c 

45 

0.51511093,-02 

-0.15172823,-02 

.62587 

.68141 

c 

46 

0.28890434,-02 

0.55486970,-02 

.74942 

.68765 

c 

47 

-0.48390226,-02 

0.15556802,-02 

.62551 

.68747 

c 

48 

-0.26852842,-02 

-0.51537342,-02 

.73422 

.67987 

c 

49 

0.46326793,-02 

-0.13908531,-02 

.62707 

.68065 

c 

50 

0.25724987,-02 

0.49868710,-02 

.74606 

.68657 
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S.  Convergence  of  the  derived  expansions 


The  estimates  obtained  in  section  4  for  rates  of  convergence  of  the 

{e  }  throw  light  on  the  convergence  of  the  series 
n 

£Cn*n<1)(2)’  Ecn*n<2><r>  <5-D 

for  the  same  four  types  of  data  : 

(i)  Smooth  data  satisfying  the  compatability  conditions 
We  recall  from  Appendix  B  that  the  norms 

II  II*  II  ^n^  II  ar*  °f  ordern^  (ln  n^3^2  (5.2) 

while  the  L  norms  are 
00 

•upyl*nl)(y)l*  l*n2)(y)l  *  °fc/ln»)  (5.3) 

^  i 

In  this  case,  by  (4.22) ,  lcnl'''  n  »  so  the  series  (5.1)  are  absolutely 

and  uniformly  convergent.  The  fact  that  they  converge  to  the  prescribed  data 

functions  f^\  f^  is  proved  by  the  author  in  another  paper  [  ]. 

Gregory  (1980)  has  given  a  proof  of  completeness  of  the  functions  <f/2^} 

n  n 

for  twice  differentiable  f^\  f^2\  with  f^  ,  f^  of  bounded  variation 
on  [0,1] . 

(ii)  Data  viol ating  compatability 

In  this  case,  from  (4.26),  |cn|  ■  0(i) ,  and  it  is  not  possible  to  guarantee 

that  the  series  (5.1)  are  convergent.  It  is  in  any  case  impossible  for 
r  (1) 

l  cn$n  (y)  to  converge  at  y  *  1  to  a  sum  different  from  zero,  since  each 
^n^(D  *  whereas  for  the  distribution  examined  in  section  4(ii),  f^(l)  ■  1, 


and  convergence  of  the  aeries  near  y  •  1  could  not  be  expected.  In  fact  the 
partial  suns  were  found  to  be  oscillatory  for  fixed  y,  but  the  work  of 
Joseph  &  Sturges  (1978)  suggests  that  the  series  should  be  Cesaro  summable  to 
the  data  as  all  points  except  y  ■  1.  This  was  confirmed  in  calculations  by 
Mayes  (1982)  for  the  corresponding  problem  for  a  cylinder. 

However  the  individual  terms  of  the  series  can  be  integrated  to  give 


0*'(t)  dt  -  (<|<  (y)  -  <M0)  ) 

n  n  n  n 


(5.4) 


s\ 


(2)(t)  -  -  <l»n1>(y)/x 

n  n 


(5.5) 


The  norms  of  these  integrals  are  asymptotically  of  order  —  times  the 
corresponding  norms  (5.2),  (5.3). 


Therefore  for  c  of  order  —  ,  the  series 
1  n1  n 


l  c«f  Jt> 
0 


dt 


(5.6) 


are  certainly  summable  in  norm,  and  computations  in  Table  7a  show  that 
with  good  accuracy  these  series  sum  pointwise  to 


fVi;(t)dt  (  -  i  y‘), 


,y  (2) 

fU'(t)dt  (-  0) 


(5.7) 


respectively.  The  quantity  a  quoted  in  the  table  is  the  least  squares  error  given 


in  each  case  by 


°2-  i 


;  -c  r 


f(t)dt  -  S(y.) 


(5.8) 


where  S  is  the  relevant  sum  fe.6  )  evaluated  at  y.  »  i  (k  -  10). 


are  absolutely  and  uniformly  sumnable.  For  the  case  computed,  with  c  ■  J,  the 
calculations  summerised  in  table  7b  show  that  the  sums  agree  with  the 

_3 

twice-integrated  data  with  error  of  order  10 


(21  * 

(iv)  Discontinuity  in  f  ' 


Similar  remarks  apply  in  this  case,  but  since  | c^ |  ■  0  In 


,  it  is 


only  necessary  to  integrate  once.  The  case  c  »  J  gives  rise  to  the  figures  in 
table  7c  showing  comparable  agreement  with  the  once-integrated  data. 
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Table  7 


Convergence  of  series  of  integrated  term* 


7a:  Data  violating  compatability  condition 


,«> .  y.  f'2> 

-  0  c  ^ 

n 

y 

j>  -  h' 

i R* 

0 

0 

0 

.1 

.005 

.00377 

.2 

.020 

.02004 

.3 

.045 

.04370 

.4 

.080 

.08016 

.5 

.125 

.12356 

.6 

.180 

.18035 

.7 

.245 

.24340 

.8 

.320 

.32055 

.9 

.405 

.40336 

1.0 

.500 

.49934 

a  -.00102 

7b: 

Discontinuity  in 

function 

f<»-  o.  f<2> 

■{-!  i«.i' 

y 

If'm 

kM2> 

0 

0 

0 

.1 

.005 

.00510 

.2 

.02 

.01999 

.3 

.045 

.04511 

.4 

.08 

.07997 

.5 

.125 

.12514 

.6 

.170 

.16995 

.7 

.205 

.20517 

.8 

.230 

.22993 

.9 

.245 

.24515 

1.0 

.250 

.25007 

(  o  *  .00010  ) 


(  J denotes  J  dy  in  all  cases  ) 


0 

-.00008 

.00016 

-.00027 

.00038 

-.00057 

.00078 

-.00112 

.00153 

-.00218 

0 

a  -  .00093 


-1/4 


0 

.00037 
.00076 
.00111 
.00153 
.00186 
.00229 
.00260 
.00304 
.00336 
.00377 

(  0  =  .00223  ) 

i 
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7c 


Discontinuity  in  first  derivative 


f(D 


-  0 


f<2>  -  1  -  8  [y-JJ 


|c  Ivi 

1  n1 


-5/4 


y 

f  (2) 

Ic  0(2) 

L  nTn 

f(2) 

Ic  *<2> 
L  n  Tn 

0 

1 

-4.463 

0 

0 

.i 

1 

-4.726 

.1 

.09994 

.2 

1 

-5.464 

.2 

.20013 

.3 

1 

-6.807 

.3 

.29980 

.4 

1 

-8.310 

.4 

.40027 

.5 

1 

-5.999 

.5 

.49963 

.6 

.2 

-.996 

.56 

.56046 

.7  ' 

-.6 

-1.023 

.54 

.53947 

.8 

-1.4 

-2.638 

.44 

.44055 

.9 

-2.2 

-5.283 

.26 

.25956 

1.0 

-3 

-3.742 

0 

0 

differs 
from  zero 
with 

O-  .00062 


(a  -  .00034) 
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Appendix  A 


1.  Stresses  and  displacements  in  terms  of  Airy  stress  function 
With  standard  suffix  notation  the  stresses  are 


“ll  "  ^ ,22 *  °12  *  -\l2- 


a22  ’  ’.11 


<A1) 


The  displacement  gradients  in  plane  strain  are  given  by 


*'*1.1  ■  <1-',)0U  -  ^22.  2<m2.2  -  •'”11*  (1-y)°22 *  ) 

“‘"l.l  ♦  “l.!5  ■  °12 


L  CA2) 


To  treat  boundary  value  problems  it  is  convenient  to  introduce  P  *  A¥, 
together  with  its  harmonic  conjugate  Q  linked  through  the  Cauchy-Riemann 
equations  P  .  ■  Q  _,  P  -  -Q  .  and  defined  so  Q  (0,0)  -  0. 

I A  9* 

The  direct  strains  are  then 


2V  “1,1  -  *,22  '  “  P’  2V  “2.2  "  *,11  ’  VP  *  "*,22+  (1’“>P 


while  the  shear  strains  are 

2“  “1.2  '  -  ’.12  *  a‘V)<> 


2“  “2.1  ■  *  ’.12  * 


Thus  if  boundary  values  are  defined  on  x^  «  0  as 


f(l) 

< 

’.12 

f<2> 

*,22 

f(3) 

■ 

Q 

,P  . 

x^  »  0 
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(A3) 


(A4) 


(A5) 


the  boundary  values  of  the  tractions  and  of  the  displacement  gradients 
with  respect  to  X2  are  respectively 


'12 


-  fd> 


:(2) 


'11 


(A6) 


3 

*  * 

U1 

-  f(1)  -  (l-v)f(3)' 

2U  ^ 

U2. 

-  f(2)  ♦  (l-v)f(4) 

.  4 

Therefore 

(i)  if  Uj  (0,x2)  and  (0,x2)  are  prescribed,  f^  and  f^ 
are  known 

(ii)  if  u2(0,x2)  and  0^(0, x2)  are  prescribed,  f^  and  f^  are  known. 

(iii)  (the  stress  problem)  :  if  o12(0,x2)  and  0^(0, *2)  are  prescribed, 

and  are  known 

(iv)  (the  displacement  problem)  :  if  Uj(0,  x2)  and  u2  (0,x2)  are 
prescribed,  the  combinations 

(A8) 

are  known. 

(i)  and  (ii)  are  canonical  problems,  (iii)  and  (iv)  are  non-canonical. . 

2.  Steady  Stokes  flow 

Neglecting  inertia  terms,  the  equation  of  motion  is  taken  as 

Vp  -  yAg  (A9) 

where  p  is  pressure  and  g  *  (u^  u2)is  the  velocity,  expressible  in  terms  of 


g(1)  ■  £«>  .  (1-v)  f<3> 


and 


g(2)  .  f<2>  -  (l-v)  £<‘> 
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the  stream  function  as  f  u2  *  1  ‘ 

Then  with  P,Q  defined  as  in  the  previous  subsection,  (  A9  )  may  be  integrated 
to 

p  »  -  yQ  +  constant 

The  stresses  are 


o. .  •  y(u.  .  ♦  u.  .)  -  p  i.. 
ij  i.J  J.J- 


whence 


2  y 


'11 


'12 


»12  »  JQ  -  p0/2y 

'*',22  '  ‘P 


Hence  corresponding  to  the  boundary  value  problems  for  the  elastic  strip, 
those  for  Stokes  flow  in  a  semi-infinite  trench  are 


(i) 

u2  and  Ou 

prescribed  on  x^ 

-  0 

—  fd) 

1  9 

f(3) 

given 

(ii) 

Uj  and  o12 

prescribed  on  x^ 

-  0 

-  £«>. 

f(4) 

given 

(iii) 

u  ^  and  u2 

prescribed  " 

.  II 

-  f«>. 

f(2) 

given 

(iv) 

°nar,d  °i2 

prescribed  '* 

•I 

-«(U. 

*(2) 

given. 

where  g(1)  -  f(1)  +  J  f(3),  g(2)  -  f(2)  -  |  f 


(4) 


(A10) 


(All) 


(A12) 


This  is  the  same  combination  of  data  as  arises  in  the  case  of  the  elastic 
displacement  problem  (A8)  with  v  ■  J. 


Appendix  B 


Quadratures  and  Norms  of  eigenfunctions 


The  quadratures  quoted  in  section  3  are  derived  from  the  following 

results,  all  of  which  are  obtained  by  integration  by  parts  and  use  of  the 

2  2  2 

equation  and  boundary  conditions  satisfied  by  the  ^  [i.e.  (D  +  X£)  ^  ■  0, 
*  D^(l)  ■  0,  D^(0)  *  D  4*^(0)  *  0, together  with  the  characteristic 

equation  2Xfc  ♦  sin  2A^  »  0] . 


1  Quadratures 

(a)  a  +  n  All  results  expressible  in  terms  of 


C  *  (cosX  coaX  )  1  /cosX  y  cosX  y  )  »  (X  tanX  -  X  tanX_)/A2  -  X2 ) 
mn  m  n  \  a'  *  /  a  a  n  n /  n ' 

He  find  (more  details  are  given  in  Appendix  A  of  [1])  : 

$(1)\  -  -AX  X  (X2  +  X2  )  C  /(X2  -  X2  )2 

Vn  Tn  /  an  a  n  an  /  a  n 

<^a2>  *n2>^  "  "  8Am  Xn  Cmn  /  <Xm  "  Xn  ^ 

\a  n  f  an  an  /  a  n 

/*(1)  <J>(3)\  -  -  AX  X  C  /  (X2 

\^m  n  /  a  n  an  •  m 

<♦”>  ♦.o,> 


- 


-A  X  C 

n  an 


/  (X2  -  X2) 

/a  n 


AC  -  A  (tanX  ♦  tanX  )  /  (X  ♦  X  ) 
an  a  nan 


(b)  a  -  n 


<♦“>  c)  *  i  •  (*?>  »:2>) 

♦  <»>  -  W  <£>♦">)  -  0 

<♦:”  ♦”’>  ■  <?i2>  c>  -  *• 


X 

m 

cosX 


nr 


The  biorthogonality  relations  (2.9)  and  (2.10)  follow  at  once  from  these  results. 

-A6- 


2. 


Noras 


The  norms  are  found  from  the  preceding  expressions  by  writing  n  ■  -m, 

so  X  «X  etc.  The  details  are  given  in  appendix  B  of  [1], 

n  m,  n  n 

and  only  the  asymptotic  expressions  for  large  m  are  quoted  here  (in  fact  they 
are  close  approximations  for  m  >  2) . 

|U  ||  *  (2am)"1  (iln4mtr)"3/2 
in 

II  II  ,  |U*2)||  *  *»"/<!  ln4aTr)3/2 

IU*3)|I»  II  II  *  2/  <ln4air)1/2 


Appendix  C  SoluCion  of  Che  biharmonic  equation  A2¥  -  0  in  the 
quarter-plane  0<8<  ~  ,  with  boundary  conditions 


e-Os^  -0,  ¥  -0 

xy  xx 


9  "  I  :  ¥xy  "  (1"V)Q  "  8<1>(y)*  Vyy  ‘  (1_V)P  "  8(2)<y> 


(g^,  g^  prescribed). 


The  biharmonic  equation  is  automatically  satisfied  when  ¥ 
has  the  form  of  a  Mellin  integral 


*  ‘ST  J*1"’ 


with  ¥  -  Acos(s-l)6  +  Bsin(s-1)8  +  Ccos(s+1)0  +  Dsin(s+1)0 


From  (C3  )  we  have 


* "  ■  m  f’1 


-1-sPds 


with  conjugate 


where 


-  -I,  f  t*l‘ 

2iri  J 


(I)  •  ■-<  (: 


(s+l)0  sin(s+l)0\  ( C 


sin(s+l)0  cos(s+l)0y  \D 


The  boundary  conditions  P  1)  are  applied  in  the  form  ¥  -0,  t* 


0-0,  giving 


A  +  C  -  0 


(s-l)B*(s+l)  D  -  0 


The  derivatives  on  8  *  y  that  occur  in  the  boundary  condition  C C2  )» 
namely 


(C8) 


are  also  expressible  as  Mellin  integrals  derived  from  (  C3  )•  When  these 
results  are  combined,  and  A  and  B  eliminated  using  (  C7  ),  the  boundary 
data  on  9  •  ^  are  expressible  in  terms  of  the  remaining  constants  C  and  D 

in  the  form 


■  &  JV‘"  s  W)-. 


with 


M(s) 


where  M  is  the  matrix 


(s-2(l-v))  cos  ^ 


(s+l-2v)  sin  ~ 


a  -  1,  2 


(s-l+2v)  sin  Tp-  ^ 
-  (s+2(l-v))  cos  —j 


(C9) 


(CIO) 


The  inverse  matrix  M  *  gives  C  and  in  terms  of  the  g  ' s,  and  the  values 
*  *  7T 

of  P,  Q  on  0  •  j  are  therefore 


(Cil) 
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(C12) 


where  A  =  -  det  M  *  s^  ♦  (3-4v)  (sin2  y  )  -  4  (1-v)^ 

The  solution  quoted  in  section  4  of  the  present  paper  is  obtained  by 
setting  v  ■  1  in  these  results.  Then 


_(1)  m  y  |  (2)  m  ^  j 

8  xy*x-0  *  8  yy*x  «  0, 


(C13) 


and 


.  2  .2  its 

A  •  8  -  Sin  y  . 


(C14) 
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Appendix  D  : 


Asymptotic  expressions  for  the  integrals 


***  dt 


(Dl) 


for  0  <  Re  v  <  1,  A  E  A  +. 

Here  A+  is  the  set  of  zeros  of  2A  +  sin  2A  lying  in  the  right  hand  half  plane 
so  that,  by  (1.10),  arg  A  ■  0  • 

In  each  case  the  integrals  are  evaluated  as  the  difference  between  two 
infinite  integrals 


...  v—1  iAt  .  V”  V* 

(l)  ft  e  dt  •  f  -  f 
0  0  1 


These  integrals  both  exist  since  Im  A  >  0.  The  first  is 
A  v  e  while  integration  by  parts  gives 


iA 

e 

1” 


C1  *  °  (x)j 


/v—1  iAt . 

'  e  dt 

Since  (as  e*X  »  o  (A  ^),  the  last  term  is  of  order  A  ^  and  altogether 


,  *v-l  iAt ...  ,-v  iviT/2r. 

ft  e  dt  ■  A  e  r(V) 


(ii)  Likewise  /  tV  *e  ^^Cdt  » 


iA 


-5/2n 


♦  ^  *  0(A  ) 


(D2) 


-l«o 

/  *  / 


and  in  the  same  manner  we  obtain 


1  v-1  -iAt, .  ,-v  -  i\rn/2 

f  t  e  dt  •  A  e 


r(v)  - 


e'iX  ’3/2 

V  +  0  (X  > 


(D3) 


The  results  quoted  in  ^.20)  are  obtained  by  addition  and  subtraction  of 
(  D2  )  and  (  D3  ) . 
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Appendix E.;  Derivation  of  equation  (4.27) 

A  direct  derivation  is  provided  by  writing  the  solution  of  the  bihartnonic 
equation  m  the  half  plane  x  ^  0  as  a  Fourier  integral 


~k-  /  *'i&  ?  (x. 


C)dC 


where 


(-)  -  c2 

Vdx/  * 


4*  -  0 


The  solution  bounded  as  x  +  “  is  f  *  (A+Bx)e~X^ 


Then  on  x  *  0, 


“hlU  'V°'t>  ’  h/  ‘C(|C|A-B)e'i5td5 


>00 

Now  the  Cauchy  integral  —  J  e~^^tdt  ■  -  i  e_i^y 

J-co  t-y 


(sgnO 


Hence  on  reversing  the  order  of  integration. 


2  f  V 


(0,t)dt 


t  -  y 


£  f  C(sgnC)  (|5|A-B)e-1^yd5 

V-00 


-  V  «>■*> 
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